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Abstract In this article we consider the static spherically symmetric spacetime 
metric of embedding class one. Specifically three new electromagnetic mass models 
are derived where the solutions are entirely dependent on the electromagnetic field, 
such that the physical parameters, like density, pressure etc. do vanish for the 
vanishing charge. We have analyzed schematically all these three sets of solutions 
related to electromagnetic mass models by plotting graphs and shown that they can 
pass through all the physical tests performed by us. To validate these special type 
of solutions related to electromagnetic mass models a comparison has been done 
with that of compact stars and shown exclusively the feasibility of the models. 

Keywords General Relativity; equation of state; electromagnetic mass; compact 
stars 


1 Introduction 

It is a widely accepted concept that the n dimensional manifold Vn can be em¬ 
bedded in a pseudo-Euclidean space of m = n(n-|- l)/2 dimensions. The minimum 
extra dimensions, m — n = n(n — l)/2 of the pseudo-Euclidean space needed is 
called the embedding class of Vn- In case of the 4 dimensional relativistic space- 
time, the embedding class is obviously 6. The well-known cosmological metric 
of EriedmannLemaitreRobertsonWalker (FLRW) [1] is of class 1, whereas the 
Schwarzschild’s interior and exterior solutions are of class I and 2 respectively. 
The Kerr spacetime metric has been shown to be of class 5 [2]. However, in the 
present paper we are limiting ourselves to the static spherically symmetric metric 
of embedding class 1 spacetime. 
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It is seen that the above mentioned metric is compatible only with two perfect 
fluid distributions, viz. (i) Schwarzchild’s solution [3] and (ii) Kohler and Chao [1] 
solution. We would like to exploit this metric to construct electromagnetic mass 
models under the Einstein-Maxwell framework by considering a charged perfect 
fluid distribution. In general, when the charge is zero in a charged distribution 
of matter, the subsequent distribution becomes the neutral counterpart of the 
charged distribution. This neutral counter part may belongs to the type of either 
Schwarzschild interior solution [3] or Kohler-Chao interior solution [3] . 

However, every charged fluid distribution indeed does not possess its neutral 
counter part and consequently if the charge is set to zero then the describing metric 
turns out to be flat and the corresponding energy density and fluid pressure will 
vanish identically. This special type of charged fluid distribution is said to provide 
an electromagnetic mass model. In connection with his model for extended electron 
Lorentz [5] conjectured that there is no other, no true or material mass, and thus 
proposed electromagnetic masses of the electron. Later on Wheeler and Wilczek 
[7] pointed out that electron has a mass without mass. Feynman, Leighton and 
Sands [8] actually termed this type of models as electromagnetic mass models. 
For further reading on historical notes and technical works on the electromagnetic 
mass one may look at the Ref. [9] and |T0l[TTl[T2l[T3l[Til[T5l[T6ll 1 711 r^n9| respectively 
under the framework of Finstein-Maxwell theory. 

Unfortunately, the electromagnetic mass models proposed by most of the above 
investigators [iiiiiiiniiiiiiaiiaiiTiiiHiDra] suffer from a negative pressure or den¬ 
sity of the fluid due to the equation of state of the form 

p + p = 0, (1) 

where p is the density and p is the pressure. This type of equation of state in the 
literature known as a false vacuum or degenerate vacuum or p-vacuum Enumm 
I23| . It has been argued that though, in general, this equation of state leads to 
negative pressure but provides easier junction conditions and realistic expression 
for mass [HIIlITlIll. Although the junction conditions do not require the density 
to vanish at the boundary as is true for gaseous spheres. Such a model is available 
in the literature for both uncharged and charged cases mm- However, we also 
note that the classical models of electron should contain the regions of negative 
density |27l[28| . It would be interesting to mention that a Weyl-type character of 
the field has been attributed which form electromagnetic mass model [5^ . 

In the present study we have attempted to obtain a charged fluid of class 1 by 
choosing specific metric potential(s) of the class 1 such that they do not form a sub¬ 
set of the metric potentials of Schwarzschild’s interior metric (inclusive de-Sitter 
and Einstein universe) and Kohler-Chao metric [3]. We argue here that the static 
spherically symmetric metric of embedding class one is more suitable to construct 
electromagnetic mass model as it possess lesser number of neutral counterparts of 
the charged fluids in comparison to general static spherical symmetric metric. Now 
if the charge be zero in the charged fluid, the describing metric will turn into flat by 
virtue of the structure of the metric. In the past, several alternatives were used by 
several investigators to obtain the electromagnetic mass models [niiiiiiiisiiii 
by employing the equation of state o as pure charge condition [13] which takes 
the equivalent form as guga = —1 m- On the other hand, Ponce de Leon [T6| 
has utilized the charged Einstein’s clusters EniEi] to get the electromagnetic mass 
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models. For further studies on different aspects of electromagnetic mass models 
one may look at the Refs. |321l33l[3Hl351l361l^l38| . 

However, for the construction of electromagnetic mass models we invoke a 
different method by adopting an algorithm which is very efficient to generate solu¬ 
tions of the desired form and physics, as such no ad hoc assumptions are required 
to obtain electromagnetic mass models. The main motivation of the present paper, 
therefore, is to obtain a set of solutions for the electromagnetic mass model with 
the help of charged fluid distribution of spherically symmetric class one metric. 
The logic behind considering the class one metric is that if one removes charge 
from the solutions then either the Schwarzschild solution [3] or Kohler-Chao solu¬ 
tion [3] will emerge the metric being flat and all the physical parameters - pressure, 
density etc. - become zero. To this aim our scheme of investigations are as follows: 
in Sect. 2 we provide the class one metric and fit the metric potentials in to the 
Einstein-Maxwell field equations for the spherically symmetric matter distribution. 
The next part is to construct electromagnetic mass models for stellar systems we 
provide the necessary algorithm (Sect. 3) and by exploiting the mathematical 
formalism we generate three new set of solutions in connection with electromag¬ 
netic mass models (Sect. 4). We also discuss the boundary conditions regarding 
all these solutions and determine the unknown constants of integration (Sect. 6). 
As till now we don’t know exact nature of the solutions set so we adopt in Sect. 
7 some specific techniques to explore the different features and properties of the 
electromagnetic mass models for physical acceptability of the anisotropic stellar 
models. In Sect. 8 we try to validate the solutions set related to electromagnetic 
mass models with some of the observed compact star candidates. We discuss our 
results in the concluding Sect. 9. 


2 The class one metric and the Einstein-Maxwell field equations 

Let us consider the static spherical symmetric metric to be 

ds^ = —e^dr^ — r^[d6^ -\- sin^Odcp^) + di^, ( 2 ) 

which may represent the space time of emending class 1, if it satisfies the Kar- 
marker condition 


i ?1414 = 


R1212R3434 + R1224R1334 
R2323 


( 3 ) 


with i ?2323 0 |30| . 

The above condition along with (ED yields the following differential equation 


Xv' 


(1 -e^) 


o " I ' 2 \ I /2 I / 

= -2[y -\-u ) Xi' +Aiz, 


( 4 ) 


with the constraint 1, where A and v are metric potentials of the line element 
ED which are function of radial coordinate r only. 

The solution of the above differential equation (BD can be obtained as 


= 1 -h A- 


( 5 ) 
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where K is non zero arbitrary constant, v'{r) ^ 0, = 1 and i'^(O) = 0. 

If Eqs. ([2]) and ([51) describe a charge perfect fluid distribution then the functions 
A{r) and v(r) must satisfy the Einstein-Maxwell field equations 

G^3 = R'j - \R9 3 = <T'3 + E' 3 )^ ( 6 ) 

where k = Stt is the Einstein constant with G = c = 1 in the relativistic geometrized 
units. 

The matter within the star is assumed to be locally a perfect fluid and con¬ 
sequently Tj and E®, the energy-momentum tensors for the fluid distribution and 
the electromagnetic field tensors, are respectively defined by 

+ p)'v''Vj - (7) 


E^3 = + Y.F'^^Fmn), (8) 

where f® is the four-velocity as = <5*4, 0® is the unit space-like vector in 

the direction of radial vector, P® = p is the matter-energy density and 

p is the fluid pressure. 

The above anti-symmetric electromagnetic field tensor Fij in Eq. (ED, denotes 
the velocity and can be defined 


®^ dAi dAj ’ 

This should satisfy the Maxwell equations 


(9) 


Rik,j “f Rkj^i F Fj^ j^ 0, 


( 10 ) 


and 


dx^ 


{V^gF ) = -A-Ky/^gF, 


where g is the determinant of quantities gij in Eq. dill) and is given by 


g = 


/e®® 0 0 

0 -e^ 0 

0 0 -r^ 

VO 




= -eG+FFsin^e, 


0 0 —r^sin^B j 


( 11 ) 


where Aj = (</)(r), 0,0, 0) is the four-potential and J® is the four-current vector 
defined by 


ji G dx^ 


where a is the charged density. 

For static matter distribution the only non-zero component of the four-current 
is and because of the spherical symmetry this has only a functional relation with 
the radial coordinate r. The only non-vanishing component of the electromagnetic 
field tensor {F^^ = —F^‘^) describes the radial component of the electric field. 
Hence, from Eq. (ED, one can easily get the expression for the electric field 


^41 _ ^ —(i/'+A)/2 



( 12 ) 
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where q{r) represents the electric charge contained within the sphere of radius r 
and is dehned by 

q{r) = Att f ar^dr = . (13) 

Jo 

Equation (US can be treated as the relativistic version of Gauss’s law which 
reduces to the following form: 

(14) 

or 

For the spherically symmetric metric (ED, the Einstein-Maxwell held equations 
can be expressed in the following ordinary differential equations 


-kT^ 


1 



r2 



(15) 


- kT'^2 = -kT^3 = 


v" \'v' v' v' - 

~2 T 


-X , <? 

e = Kp + — 
^4 


^4 _ A' _A , (1-e-^) _ , q^ 

— k1 4 — —e H- 5 - — Kp —j, 

r ^ r^ 


(16) 


(17) 


where the prime denotes differentiation with respect to the radial coordinate r. 
By using the Eqs. (USD, (USD, e 3 and also we obtain 


- - - -t; -r(4;r — Ku') = Kp — J-r, 

{A +Ku'^e’^y 


(18) 


(4 + Kiy'^e'') 

Kv'e'^ 


p' [Kv'e'' - 2r) {2v" + i/'^) 
2r “ ' 


2r (4 + Kv''^e'') 


= Kp+-r, 


(4 + Ku'^e’') 


A{2u" + u'^) ^ 
(4 + Kv''^e'^) r 


= Kp + 


On the other hand, the pressure isotropy condition can be given by 


(kv'e'^ \ 

2v' 

A{2v" + ' 

_ 2g" 

1 2, V 

r (4 -|- Kv''^e'^) 

(4-h Ki/'2e^)^ 



(19) 

( 20 ) 


( 21 ) 


A closer observation of the above set of differential equations easily indicates 
that if charge vanishes in a charged fluid of embedding class one, then survived 
neutral counterpart will only be either the Schwarzschild [3] interior solution (or 
its special cases de-sitter universe or Einstein’s universe) or the Kohler-Chao [3] 
solution, otherwise either the charge cannot be zero or the survived space-time 
metric will become hat. 

Now, one can look at Eq. EID which immediately indicates that in absence of 
the charge either of the two factors on the left hand side has to be zero. Conse¬ 
quently, it can be shown that if the hrst factor of Eq. dUD be zero then it gives 
rise to the Kohler-Chao [3] solution in the form: 

ds^ = — ^ ~ + sin^6d(fJ) -|- (A -|- Br‘^)dt^. 


( 22 ) 
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in this model, are 


(A + 2Br2) ’ 

(23) 

^{3A + 2Br^) 

= ^ (A + Br2) • 

(24) 


One can observe from (HI that since it does not possess zero pressure as well 
as density for any finite radius on the surface, it cannot represent a compact star. 

Let us now consider the second factor of Eq. dUD which in its vanishing form 
ultimately provides the Schwarzschild [H] interior solution 




— r(d6^ + sin6d(j> ) + I A B \ 1 — 




dt 


with its pressure and density as follows: 

A + ZB^Jl 


W 


up = -- 


a + bJi- 


(25) 


(26) 


Kp = 




(27) 


where A and R are non-zero constant quantities and B > 0. 

If the mass function for electrically charged fluid sphere is denoted by m(r), 
then it can be dehned in terms of the metric function as 


5-A(r) ^ ^ _ 2m(r) q 


(28) 


where the function m(r) represents the gravitational mass of the matter contained 
in a sphere of radius r. Now, if R represents the radius of the fluid sphere then it can 
be shown that m is a constant with m(r = R) = M outside the fluid distribution 
where M is the gravitational mass. Following the work of Florides |30| this can be 
defined as 

M = pi{R)+i{R), (29) 

where p.{R) = f f^pr^dr is the mass inside the sphere, ^(R) = ^ arqe^^^dr is 
the mass equivalence of the electromagnetic energy of distribution and Q = q{R) 
is the total charge inside the fluid sphere. 

By using Eq. dMl) one can write the mass, in terms of energy density and charge 
function, as follows: 


m{r) = -Jprdr + -J ^dr -f —, 


(30) 


Again from Eqs. m and (nm we obtain expression for metric potential 


' 1 - 


■ + 


( 31 ) 
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Also, the expression for the pressure, in its gradient form, can be obtained by 
using Eqs. (US and (HU) - dsni) as follows 


dp _ Ma{r){p + p) (x-v)/2 , q dq . , 

dr dr ’ 

where Mq is the gravitational mass within the sphere of radius r and is given by 

Mcir) = (33) 

The above Eq. (1321) represents the charged generalization of the Tolman-Oppenheimer- 
Volkoff (TOV) equation of continuity for perfect fluid stellar system |411I42| . 


3 Algorithm for electromagnetic mass models 


We are now in position to construct ‘electromagnetic mass’ models of stellar system 
for class one metric by using algorithm given by Maurya et al. |43| . 

The Eqs. (1151) - (1171) in terms of mass function reduce to 


2m{l-\-rv') ^ v' ^ ^ _ 

■3 I “T A I 4 ^Pl 

rpxj ly* r^ r^ 


(34) 


1^(2 + ru') m{2r^v" + r^u'^ + rv' — 2) 


+ 


2r2 2r^ 

2rqq' v' — 2q^i/' + Aqq' + (r^ + q^){2rv" + rv'^ + 2v) 2(f' 

2m' 2qq' 2 
“72-3“ = P- 


= (35) 


,0 (36) 

From Eqs. dMI and (ESI), the first order linear differential equation for m(r) in 
terms of v{r) and electric charge function q{r) can be provided as follows: 


, (2r^v" -\-v'^ — “iu'r — Q) r(2rv"rv''^ — 2v') , . 

m ^^- -m = — -^77—r—777-^ + /(»') 


r[ru' + 2) 


2{rv' + 2) 


where 


^ q^ [2r'^v" + rv'{rv' - 4) - 16] qq'{rv' + 2) 

2r'^(rv'-\-2) r{rv'-\-2) 

Hence the mass function m(r) can be given by 

m{r) = e J {/i(r) +/(r)} dr + A 

{2r^v" + r^v''^ — 3rv' — 6) 


where 


g{r) = 


h{r) = 


r{rv' + 2) 

r{2r^v" + rv''^ — 2v' 
2{rv' + 2) ■ 


(37) 

(38) 

(39) 

(40) 


and 


(41) 



















4 New class of electromagnetic mass models for stellar systems 


4.1 Solution of Type I 

We consider the following suitable function 

= 2Ar^ + logB, (42) 


A(r) = log 1 + K 




where A and B are positive constants. 

The expressions for the mass and the electric charge are respectively 


(43) 




De 


2Ar-‘ 


1 + DAr‘^e‘^^^‘ 


+ 


+ A - ) 

2(l + DAr2e2^’-")2 


(44) 


^=E^= A^r^ 

r'4 


4 _ 475e2^’'" 
2(l + DTr2e2^’'")2 


(45) 


where D = AABK is a pure constant. 

Again, the expression for the energy density and the pressure are given by 


Svrp = A 


D^Ar^e^^^ - 4Ar^ + 6De^^’' (2Ar^ + 1) 

2(l + DAr2e2^’-")2 


(46) 


Sttp = A 


-D^Ar‘^e'^^^~ + 4(2 + Ar^) + 2De^^'"''{2Ar^ - 1) 


2Ar'’ 


2(l + DAr2e2^’-")2 

The respective gradients of above physical parameters are 


(47) 


dp 

dr 


2A^r -D^Ar^e^^^ + (-3 + AAr^ + SA^r^) + (4 + 7Ar^ + dAV^) - 4 


S-TT 


2(l + DAr2e2^'-")3 


(48) 


^ _ 2A^r 
dr 8tt 


4)3Ar2g6Ar" D2g4Ar"(4^ 20Ar2 + 28A2r^) - 4De2^’’"(6 + TAr^ + 4AV^) + 4 

2(l + DAr2e2^’-")3 


( 49 ) 
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4.2 Solution of Type II 

Here the functional relation for the metric potentials are 

v{r) = 2log{l + s\nh Ar^) + logB, (50) 


= [l + K 




where A and B are positive constants. 

The expressions of mass and electric charge are 


(51) 


MO =^ 1^2 


D cosh^ Ar^ 


—Ar 


r (1 + cosh^ Hr^) 

2 f D{2 sinhTr^ coshHr^(l + sinhHr^) + 2 cosh^ Ar^) — cosh^ Ar^(l + sinhHr^) — 4sinhTr^ 


2(1 + DAr'^ cosh^ Ar‘^)'^{l + sinhHr^) 


(p2) 


^2 o o / —H r(2sinhTr^ coshTr^ — Z)cosh^ylr^)(l + sinhHr^) + 2cosh^ +4sinhTr^ \ 

S_ = E^ = A^r^ -^- Al - i -^-(ks) 

y 2(1 + cosh^ Hr2)2(l-|-sinhHr^) J 

The energy density and pressure (taking x = Cr^) are given by 

—Ax sinh x + D [2a; cosh^ ® + (6 cosh^ x + 10® sinh x cosh x + Dx cosh"^ ®) (1 + sinh ®)] 


Snp = A 


2(1 + Dx cosh^ x)2(l + sinhx) 


(54) 


® cosh^ ® — (2 cosh^ X + 2® sinh® cosh® + D® cosh^ ®)(1 + sinh® 
2(1 + D®cosh^ ®)^(1 + sinh®) 

(55) 

As done in the previous case, the expressions for the pressure and the density 
gradient can be determined by taking their derivatives with respect to r which are 
not produced here being their very complicated forms. 


Sarp = A 


8 cosh ® + 4® sinh ® + D [6 


4.3 Solution of Type III 

The metric potentials in this case are related to the following functions 

ia(r) = 2log(l + sin Ar^) + logB, (56) 


gA(r) ^ 


1 +A- 


(57) 


where A and B are positive constant. 



10 



r/R 


Fig. 1 are plotted with continuous line for solution I, small dashed line for solution II and 
long dashed line for solution III. 



r/R 


Fig. 2 are plotted with continuous line for solution I, small dashed line for solution II and 
long dashed line for solution III. 


The expressions of mass and electric charge are 


_ 4,,2 

r 


D cos^ Ar^ 


(1 + DAr'^ cos^ Ar'^) 

. 2 f -D(sin2j4r^(l + sin — 2 cos^ Ar'^) + cos^ Ar'^{l + sin Ar^) — 4sin \ , 
^ ^ 2(1 + cos^ Ar2)2(l-|-sin^r^) y 



.2 2 

= A r 


' D [sin 2Ar^ (1 + sin Ar^) — 2 cos^ j4rcos"^ Ar^(l + sin j4r^) — 4 sin Ar'^ 


2(1 + DAr'^ cos^ Ar‘^)‘^{l + sin^r^) 
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The expression for energy density and pressure (taking x = Cr^) are respec¬ 
tively 


Srrp = A 


4* sin a: 4- D [23:cos^ a; -|- (6 cos^ x — 5a; sin a;) (1 + sin a;)] -|- D^a:cos"^ a:(l -|- sin a;) 


2(1 + Dx cos2 a:)2(l -|- sina;) 


(60) 


8 cos a; — 4a; sin x 
Sttp = A - 


+ D [ 6 a; cos^ a; — (2 cos^ x — x sin 2x 4 - Dx cos^ *) (1 + sin a;)] 
2(1 + Dx cos^ a;)2(l 4 - sina;) 


(61) 

Likewise the case II, the expressions for the pressure and the density gradients 
being very cumbersome we are leaving those calculations of respective derivatives. 

Let us look at Figs. 1 and 2 regarding the desirable features on the basis of their 
respective solution. It is expected that the solution should be free from physical 
and geometrical singularities, i.e. the fluid pressure and the energy density at the 
center should be finite and metric potentials and should have non-zero 
positive values in the range 0 < r < i?. At the center one must have = 1 
and = B for each solution. Interestingly, both Figs. 1 and 2 show that metric 
potentials are positive and finite at the center. 

Similarly, the density p should be positive and the pressure p must be positive 
inside the star as well as it should be zero at the boundary of the fluid sphere. All 
these features are quite available from Figs. 3 and 4. 



r/R 


Fig. 3 Pressure is plotted with long dashed line for solution-I, continuous line for solution II 
and small dashed line for solution III. 


Let US summarize the above results and at a glance try to get the flavor of these. 
We would like to mention here that if A = 0 in the cases (I), (II) and (III) then 
the corresponding metrics at once turn to flat spacetime and also the expressions 
for the electric charge, the pressure and the energy density automatically vanish. 
Therefore, the three charged fluid distributions obtained above depict the three 
electromagnetic mass models. 
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r/R 


Fig. 4 Density is plotted with long dashed line for solution I, small dashed line for solution 
II, continuous line for solution III. 


5 Boundary conditions for the spherical system 

The above system of equations has to be solved under the condition that the 
radial pressure p = 0 at r = a (where r = a is the outer boundary of the fluid 
sphere). The interior metric ([2]) can join smoothly at the surface of spheres to the 
Reissner-Nordstrom metric [H] 

= + dr^-r^(d6»^+sin^6ld().^)+^1 - ^ (62) 


This requires the continuity of e^(r), e'^{r) and q{r) across the boundary r = R. 


e-A(fl) 


1 - 


2M Q 


R 




(63) 




2M 

R ^ R^ 


(64) 


q{R) = Q, 


(65) 


P{r=R) — 0. 


( 66 ) 


By using all the above boundary conditions we are able to find out expressions 
for various constants as can be seen below. 
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5.1 Solution of Type I 


D = 


{2AR^ - l) + i/(8l2^M^4ARM^ 

AR2^2AP 


At the boundary: 


g-A(«) = ^^(R) 


gives 


B = 


g2Afi2(i_^^^^2g2Afl2)- 


(67) 

( 68 ) 

(69) 


5.2 Solution of Type II 


The pressure is zero on the boundary r = R (here X = AR^) and we obtain 
expressions for the constants as follows: 

A(1 + sinh A)(2cosh A+ Asinh X) - 
—Acosh^ X(1 + sinh A) ’ ^ ^ 

where 

D\{X) = cosh^ A(1 — 3A cosh A + sinh A) + A cosh A sinh A(1 + sinh A) (71) 
and 


D 2 (A) = - cosh^ A (4A - 2 cosh A + 2A cos 2A - 2A sin A - sinh 2A)^ . 
At the boundary: 

g-A(fl) = 

gives 

_ 1 _ 

(1 + sinhA)^ (l ADAcosh^ A) ' 


(72) 

(73) 

(74) 


5.3 Solution of Type III 

The pressure being zero on the boundary r = R (here A = AR^) the constants 
can be given by 

„ cos^A(6AcosA-2sinA-2) + Asin2A(l + sinA) 
- 2Acos4A(I + sinA) -’ 


t ^ - I6A cos4 A (1 + sin A) (-2 cos A + A sin A) + [6A cos3 A - (2 cos^ A 

-2A cos"^ A(1 + sinA) 

At the boundary: 

g-A(«) = ^HR)^ 

gives 

B = - ^ -. 

(1 + sinA)^ (1 + DAcos2 A) 


Asin2A)(l + sinA)]2 

(76) 

(77) 


( 78 ) 





14 


6 Physical features of the electromagnetic mass models for stellar systems 

In the solution part (Sect. 4) we have analyzed some of the physical parameters, 
potentials, density, pressure etc., through their graphical plots. They exhibited 
desirable physical features regarding stellar configuration. However, in the present 
Sect. 6 we are interested to perform a few rigorous tests for other physical param¬ 
eters, velocity and charge, and also prepare a check list for energy conditions and 
stability issues (such as TOV equation and Buchdahl condition). 


6.1 Sound velocity 

The velocity of sound should monotonically decrease away from the center and 
increase with the increase of density, i.e. ^ (^) < 0 > 0 0 < r < R. 

It is argued by Canuto that the equation of state at ultra-high distribution of 
matter the sound speed decreases outwards. 

In the present model, from Fig. 5, it is clear that velocity is decreasing for 
solution I and increasing for solution II and III throughout the star. Therefore, 
the solutions for solution II and III are not suitable at all as far as compact star is 
concerned. This is because the equation of state for nuclear matter shows a regular 
behavior of ^ for these solutions |46| . 

The above discussions, based on the demonstration of the hgures, immediately 
restraint ourselves to study henceforth only the solution of type I electromagnetic 
mass model. 



r/R 


Fig. 5 Velocity is plotted with continuous line for solution I, long dashed line for solution II 
and small dashed line for solution III. 


6.2 Electric charge for solution I 

From the present model it is observed that in the unit of Coulomb, the charge on 
the boundary is 1.15295 x 10^° C and at the center it is zero (as the charge on 
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the boundary is 0.9889 so we have to multiply this by the number 1.1659 x 10^° 
to obtain the resultant numerical value). 

One can observe from Fig. 6 that the charge profile starts from a minimum and 
acquires the maximum value at the boundary. This figure has been drawn for the 
compact star RX J 1856 — 37 with the constant values CR^ = 0.1836, D = 2.9540. 
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Fig. 6 Electric Charge is plotted with long dashed line for solution I. 


6.3 Energy Conditions for solution of type I 

For physical validity an energy-momentum tensor has to obey the following energy 
conditions: 

1. null energy condition (NEC): p + ^ > 0, 

2. weak energy condition (WEC): p — p + ^ >0, 

3. strong energy condition (SEC): p — >0. 

We have plotted the feature of different energy conditions in Fig. 7 for the 
values of different physical parameters connected to energy conditions for the 
constants: CR^ = 0.1836, M = 0.9041 Mq, R = 6.006 Km and ^ = 0.222. 
The figure indicates that all the energy conditions are satisfied throughout the 
interior region of the stellar system. 


6.4 Generalized TOV equation for solution of type I 

We write the generalized Tolman-Oppenheimer-Volkoff (TOV) equation [37] in the 
following form: 


MG{p + Pr ) ^( X -^)/2 _ ^ = 0 

J.2 J.2 


( 79 ) 
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r/R 


Fig. 7 NEC is plotted with long dashed line, WEC is plotted with continuous line and SEC 
is plotted with small dashed line for solution I. 


where Mq is the effective gravitational mass within the radius r and can be pro¬ 
vided 

Mcir) = V . (80) 

The above TOV equation describes the equilibrium condition for a charged fluid 
subject to gravitational {Fg), hydrostatic (T^) and electric {Fe) forces. Therefore, 
one can write it in a more suitable form 


Fg + Ff^ + Fe — 0, 


where 


Fg = -^^'(P + P) 


2A^r (1 -f -f 4 

8n (1 -hDy4r2e2^r2)2 


(81) 


(82) 


F - ‘^P- 


and 


Stt 


+ AAr"^ -f 8TV‘‘) + + lAr"^ -h dTV'^) - 4 


2(1 + Dy4r2e2^r2)3 


(83) 


P q A/2 

= = — 


(-2 -h ) -6 -h D^Ar^e^^^ + (3 -h 2Ar^ + 8A^r^) 


2(l-hDAr2e2^r2)3 


(84) 

The plot for the TOV equation is shown in Fig. 8. We observe from this figure 
that the system under the joint balancing action of the different forces, e.g. gravi¬ 
tational, hydrostatic and electric to attain an overall static equilibrium. However, 
from Fig. 8 it is also clear that the gravitational force has the dominant role over 
the hydrostatic force whereas the electric force has the negligible contribution to 
the equilibrium. This feature seems quite reasonable in the case of the compact 
stellar system. 
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Fig. 8 Fg is plotted with long dashed line, is plotted with continuous line and Fe is plotted 
with small dashed line for solution I. 


6.5 Effective mass-radius relation and surface redshift for solution of type I 


Buchdahl |48| has proposed an absolute constraint of the maximally allowable 
mass-to-radius ratio {M/R) for isotropic fluid spheres in the form 2M/R < 8/9. 
However, Bohmer and Harko |49| have shown that for a compact object with 
charge, Q{< M), there is a lower bound for the mass-radius ratio 


3Q^ 

2i?2 


1 + _ 2 L\ 

^ ^ 18 fl^ I 

^ 4- 

^ + TTW / 




(85) 


whereas the upper bound of the mass-radius of a charged sphere was generalized 
by Andreasson |5n| as follows: 


Vm<^ + 


R 0^ 
9 3R’ 


In the present model, the effective gravitational mass is given by 


( 86 ) 


M, 


eff 



r'^dr = = ^R 


1 -h DAR'^e^^^^ 


. (87) 


Therefore, the compactness factor can be written as 

_ Me// _ 1 r DAR^e‘^^^'‘ 

R “2 TT'dAR2 ^2A^ 


( 88 ) 


The surface redshift in connection with the above compactness is given by 

Z = {l-2u)~^/‘^ -1 = - 1 = a/iT'dAR2^2Ar7_ ^ 

The plot of the surface redshift is shown in Fig. 9 for the compact star RX J 1856— 
37 with the constant values CR^ = 0.1836, D = 2.9540. It can be observed that 
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there is a gradual increase in the redshift which is an acceptable physical fea¬ 
ture. The maximum surface redshift for the present stellar configuration of radius 
R = 6.006 Km turns out to be Z = 0.3882 which seems well within the limit 
Z < 2 [1811^1521 . 



Fig. 9 Red shift is plotted with long dashed line for solution I. 


7 Validating the model with strange star candidates 

In the foregoing Sect. 6 we have studied several physical behavior of stellar system 
in connection with electromagnetic mass models. In some of the subsections, e.g. 
6.2 (electric charge) and 6.5 (surface redshift), we have also shown graphical plots 
specifically for the compact star RX J 1856 — 37 with the mass M = 0.9041 Mq 
and the radius R = 6.006 Km. 

However, it seems that some more investigations are needed to show the va¬ 
lidity of our models for other compact stars which have definite observed physical 
features. In the following two Tables 1 and 2 we, therefore, produce data sheet 
for the purpose of comparison between the present model stars and the observed 
compact stars. 


Table 1 Values of the model parameters A, B, D and K for different strange stars 


Strange star 
candidates 

M (Mq) 

R (Km) 

M/R 

B 

A 

D 

K 

RX J 1856 - 37 

0.9041 

6.006 

0.222 

0.5189 

5.0962 X 10-13 

2.9540 

2.7925 X 1012 

Her X - 1 

0.9825 

6.700 

0.216 

0.5552 

3.6319 X 10-13 

3.0626 

3.7972 X 10i2 

PSR 1937 -1- 21 

2.1 

11.4998 

0.269 

0.4103 

1.9503 X 10-13 

2.5857 

8.0775 X 10i2 

PSRJ 1614 - 2230 

1.97 

11.3664 

0.2553 

0.4419 

1.8112 X 10-13 

2.6998 

8.4338 X 1012 

PSRJ 0348 + 0432 

2.1 

11.7372 

0.2636 

0.4228 

1.8017 X 10-13 

2.6315 

8.6369 X 10i2 
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Table 2 Energy densities and pressure for different strange star candidates for the above 
parameter values of Table 5 


Strange star 
candidates 

Central Density 
{gm / cm~^) 

Surface density 

Central pressure 
(dyne ! 

RXJ 1856 - 37 

Her X -1 

PSR 1937 -1- 21 
PSRJ 1614 - 2230 
PSRJ 0348 + 0432 

2.4252 X IQl® 
1.8869 X IQi® 
8.1241 X IQi-* 
7.4524 X IQi-* 
7.6379 X IQi^ 

1.6183 X IQi® 
1.2718 X IQi® 
5.0473 X IQi-* 
4.9876 X IQi-* 
4.7797 X IQi"' 

2.2243 X 10®® 
2.5768 X 10®® 
1.3334 X 10®® 
1.1384 X 10®® 
1.1918 X 10®® 


For our model we particularly note that for the compact star RXJ 1856 — 37 
with mass M = 0.9041 Mq and radius R = 6.006 Km the surface redshift turns 
out to be .^ = 0.3882 which seems falls within the range Z < 2 HHlEIlEa and 
0 < Z < 1 [53115411551I561I57| . However, one may figure out the surface redshifts for 
other compact stars also as provided in Tables 1 and 2 and we expect those values 
will be within the above specified range. On the other hand, surface density as can 
be seen from Table 2 is of the order of 10^^ — 10^® gm/cc. This very high density 
indicates that the model under ‘electromagnetic mass’ represents an ultra-compact 
star [58115^160] . 

Therefore, we would like to pass a general remark that our models in connection 
with ‘electromagnetic mass’ represent compact stars of several categories. 


8 Conclusion 

We have considered the static spherically symmetric spacetime metric of embed¬ 
ding class one in the present investigation. Is has been possible to show the ex¬ 
istence of electromagnetic mass models specifically in connection with compact 
stars. Three new electromagnetic mass models are derived where the solutions are 
entirely originating from the electromagnetic field, such that the density and pres¬ 
sure like physical parameters do vanish for the vanishing charge alone. However, 
a meticulous analysis reveals that among these three sets of solutions all are not 
equally interesting as far as astrophysical several aspects are concerned. To val¬ 
idate these special type of solutions related to electromagnetic mass models, we 
have also conducted a comparison between our proposed model and the observed 
compact stars which shows satisfactory results in favor of the present theoretical 
modeling. 

However, an obvious question may arise to the study of the compact stellar 
configuration under Einstein-Maxwell spacetime, especially how the charge comes 
in the consideration of such kind of systems. A brief historical note on the issues of 
stability of static spherically symmetric stellar systems and as an effective measure 
of averting singularity why one should include charge and what is the process of 
holding huge amount of charge inside the bodies are available exhaustively in the 
Ref. |35| and the Refs, therein. As an continuation of this discussion, we feel, an 
outline on the charged bodies may be helpful to the readers which follows in the 
next paragraph. 

In the history of general relativity the first ever exact solutions of the Einstein 
field equations, the well-known Schwarzschild interior solutions, suffer from the 
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problem of singularity due to gravitational collapsing of a spherically symmetric 
matter distribution. One way to overcome this singularity is to include electrical 
charge to the neutral bodies. It is suggested by the scientists that gravitational 
collapse can be avoided in the presence of charge where the gravitational attrac¬ 
tion is counter balanced by the electrical repulsion in addition to the pressure 
gradient |611l621l6^ . To this end questions came up regarding the stability of the 
charged sphere and also about the amount of charge that holds by the star. A good 
amount of works have been done by several authors on the stability issue [6ill65l 
[66ll67l[68ll69[ . On the other hand, in some recent studies [ToirmiiT] we find out 
estimate of electric charge in the compact stars which amounts a huge charge of 
the order of 10^^ — 10^° Coulomb. 
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